function res = multimid(fun,xl,xr,n)
% MULTIMID - 用复合中点公式求积分。
%   input:
%   对函数fun 从 xl到xr进行积分。

    res =0;
    h = (xr-xl)/n;
    points = xl:h:xr;
    points = points-h/2;
    F = fun(points);
    res = h*sum(F(2:n+1));
end